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We propose a model of inflation based on a simple variant of the NMSSM, called 
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0NMSSM, where the additional singlet plays the role of the inflaton in hybrid (or 
inverted hybrid) type models. As in the original NMSSM, the 0NMSSM solves the \i 
problem of the MSSM via the VEV of a gauge singlet N, but unlike the NMSSM does 
not suffer from domain wall problems since the offending Z 3 symmetry is replaced 
by an approximate Peccei-Quinn symmetry which also solves the strong CP problem, 
and leads to an invisible axion with interesting cosmological consequences. The PQ 
symmetry may arise from a superstring model with an exact discrete Z 3 x Z 5 symmetry 
after compactification. The model predicts a spectral index n = 1 to one part in 10 12 . 
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1 On leave of absence from 



There is to date no standard model of inflation, and although there has been a 
good deal of progress in recent years in this area much of the current activity has 
been concerned with conceptualised field theoretic models rather than well motivated 
particle physics based models [T] . Possibly the best motivated particle physics model 
beyond the standard model is the minimal supersymmetric standard model (MSSM). 
However the only Higgs fields in the MSSM are the two doublets Hi, H 2 , which develop 
vacuum expectation values (VEVs) of order the weak scale, and it is very difficult 
if not impossible to develop a model of inflation using only these fields for several 
reasons. The primary reasons are that the electroweak scale turns out to be too 
small and the Higgs potential is not sufficiently flat. The so called next-to-minimal 
supersymmetric standard model (NMSSM) is more promising from the point of view 
of inflation since it contains, in addition to the two Higgs doublets, a Higgs singlet N 
which may develop a large VEV. 

The usual NMSSM does not require the fiHiH 2 term of the MSSM, replacing it 
with a XNH1H2 term, and thereby solving the fi problem^. The NMSSM also involves 
a term kN 3 in the superpotential so that the model has an exact Z3 symmetry ||. ff| . 
However this is broken at the weak scale leading to a serious domain wall problem 
||. U . Originally it was thought that the Z 3 may be slightly violated by Planck scale 
operators, leading to a pressure term that removes the walls. However without an 
exact Z 3 symmetry supergravity tadpole diagrams will lead to a large singlet mass in 
the low energy theory, and the amount of Z% breaking required to solve the domain 
wall problem is in conflict with requirement that tadpoles do not make the singlet 
too heavy |6|, |7|]. 

It transpires that, without fine-tuning, the NMSSM does not lead to a sufficiently 

flat potential along which the inflaton may roll. In order to overcome this we introduce 

2 Note that the Giudice-Masiero mechanism presents a solution to the \i problem within the 
MSSM by generating the [i term via a non-minimal Kahler potential. 
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a second singlet 0, and replace the term N 3 in the NMSSM by (f)N 2 . Thus our model 
is based on the superpotential: 

W^nmssm = XNH 1 H 2 - kct>N 2 (1) 

Note that our model has the same number of dimensionless couplings as the original 
NMSSM, and we have used the same notation A, k to emphasise this. With this 
modification the field <fi appears only linearly in the superpotential and so will have 
a very flat potential, lifted only by a tiny mass of order electronvolts, and will 
play the role of the inflaton field of hybrid inflation || |9|, [H| if m\ > or inverted 



hybrid inflation [n] if mi < 0. In the case of inverted hybrid inflation the present 



model provides an interesting counter example to the problems raised in Ref. [12 1 . 



Inflation ends when reaches a critical value <p c ~ 10 13 GeV after which the N field, 
which has a zero value during inflation, develops a VEV < N >~ C . Interestingly 
the inflaton also develops an eventual VEV < >~ C via a tadpole coupling, which 
is typical of inverted hybrid inflation but quite extraordinary for hybrid inflation. The 
resulting dimensionless couplings are X,k ~ 10~ 10 , whose smallness will be explained 
by embedding the model into a string inspired model where the couplings result from 
higher dimension operators, controlled by discrete symmetries. Note that radiative 
corrections to the inflaton mass are controlled by A, k and are of order the inflaton 
mass itself. 

Having replaced the NMSSM superpotential by Eq. ([!]), the troublesome Z3 sym- 
metry is replaced by a global U (1) pq Peccei-Quinn symmetry where the global charges 
of the fields satisfy: 

Qn + Q Hl + Qh 2 =0, Q + 2Q N = 0. (2) 

with the quark fields having the usual axial PQ charges. The global symmetry forbids 
additional couplings such as N 3 , <pHiH 2 and so on, but is broken at the scale of the 
VEVs releasing a very light axion. The axion scale f a is therefore of order C in this 



model. The axion will be an invisible Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) |p~3| 
type axion, which couples to ordinary matter through its mixing with the standard 
Higgses after the electroweak phase transition. Once we embed our model into a string 
motivated model, the global PQ symmetry will emerge as an approximate accidental 
symmetry of an underlying discrete symmetry, and we need to discuss such questions 
as the solution to the strong CP problem in this wider context. Note that if we had 
simply removed the N 3 term from the NMSSM superpotential and not replaced it 
with anything then the theory would also have a PQ symmetry, and the potential 
would also be flat in the N direction, and then one might be tempted to identify iV 
with the inflaton of hybrid inflation. However in such a scenario the height of the 
potential during inflation would be of order 1 TeV, leading to an inflaton mass very 
much smaller than the radiative corrections to its mass of order 1 eV, which would 
require unnatural fine-tuning. By contrast, with the cf)N 2 term present, the height 
of the potential during inflation is about 10 8 GeV and the COBE constraint may be 
satisfied by an inflaton mass of about 1 eV which is the same order as the radiative 
corrections to its mass, leading to a natural scenario with no fine-tuning required. 

The tree-level potential which follows from the superpotential in Eq. (TJ) can be 
written, if we ignore Hi, H2 which have smaller VEVs, 

V = V(0)+V(<f>,N) 
V(<j>,N) = k 2 N 4 + m 2 (<f))N 2 + ml<i)\ (3) 

with the field dependent N mass given by, 

m 2 (0) =m 2 N -2kA k (t) + Ak 2 (t) 2 . (4) 

We have taken and N to be the real components of the complex singlets, and 
included the soft breaking parameters from the soft supersymmetry breaking potential 
terms m^N 2 , m^cj) 2 and Akk<pN 2 . We have also added by hand a constant vacuum 
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energy V(0) to the potential, about which we shall say more later. Note that m 2 (0) 
for equal to a critical valueQ 




(5) 



In order to discuss inflation we need to specify the sign of the inflaton mass squared 
m\. If m 2 , > (as in hybrid inflation) then, for > 0+, iV will be driven to a local 
minimum (false vacuum) with N=0. Having a positive mass squared, will roll 
towards the origin and m 2 (0) will become negative once the field reaches After 
that, the potential develops an instability in the N=0 direction, and both singlets 
roll down towards the global minimum, 

<^> = * (6) 



<N> 



4^f = V2 c ± -<0> , (7) 



signaling the end of the inflation. On the other hand if m 2 < (corresponding to 
inverted hybrid inflation) then we shall suppose that during inflation < 0~ , with 
the inflaton rolling away from the origin, eventually reaching 0~ and ending inflation 
with the same global minimum as before. Note that the global minimum VEV < > 
is sandwiched in between <p~ and 0+ so either hybrid or inverted hybrid inflation is 
possible in this model depending on the sign of m|. 

Since Ak is a soft SUSY breaking parameter of order 1 TeV we have the order of 
magnitude results: 

fc0^ ~ k < N >~ k < >~ 1 TeV. (8) 

Since the VEVs are associated with the large axion scale, we see that the parameter 
k ~ O(10~ 10 ). Similarly since A < iV > plays the role of the \x parameter of the 
MSSM we require A to have a similarly small value. We shall discuss the origin of 



3 We require that the condition A| > ^m? N is fulfilled 
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such a small values of A, k later in the context of the string motivated model, but for 
now we simply note their smallness and continue. 

The negative value of V(4>, N) at the global minimum, is compensated by V(0) 
which is assumed to take an equal and opposite value, in accordance with the observed 
small cosmological constant. Thus we assume: 

V(0) = -V(< (j) >,< N >) = k 2 < N > 4 = 4A; 2 (0±- < >) 4 . (9) 

During inflation we may set the field N = so that the potential simplifies to: 

V = V(0) + ml4> 2 (10) 

The slow roll conditions are given by: 

The subscripts "N" means that <fi and e have to be evaluated N e-folds before the end 
of inflation, when the largest scale of cosmological interest crosses the horizon^, that 
is, N ~ 60. The height of the potential during inflation is approximately constant 
and given by 1/(0)3 = k\ < N >~ 10 8 GeV. 

Assuming that V(0) dominates the potential during inflation, <p N = (j)fe vN ps <fif, 
where the last approximation follows since in our model \rj\ <C 1/N. We need further 
to check that our inflationary model is able to produce the correct level of density 
perturbation, responsible for the large scale structure in the Universe, accordingly to 

the COBE anisotropy measurements. The spectrum of the density perturbations is 

The required number of e-folds is roug hly 60 for a potential barrier y(O) 1 / 4 ~ 10 16 GeV and 
very efficient reheating in the post-inflationary period, ft diminishes when the value of V(0), and/or 
the reheating temperature decrease. This will be the case for this model, where we will see that 
the needed value of V(0) is lower than 10 11 GeV, and the reheating temperature is quite low. 
Nevertheless, because \r)\ is extremely small, the particular value of N is not relevant to this stage. 
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given by the quantity ]14 



,2 32^(0) 1 



^ = 75M^' (13) 



with the COBE value, 5 H = 1-95 x 10 5 Hl5fl. Writing (fif ~ C , COBE gives the order 



of magnitude constraint: 

This, in turn, is more than enough to broadly satisfy the slow-roll conditions. In 
particular, 



M 2 P \k m<j> \ 2 
8tt (V2k(j) c )* 



\Vn\ - ^T^T7-10- 12 , (15) 



K \km^ A2 fa 2 

eN ~ ^(Tl^ ^ 477 ^ (16) 

The model predicts a very flat spectrum of density perturbations, as usual in this type 
of hybrid model, with no appreciable deviation of the spectral index, n = 1 + 2r] — 6e, 
from unity Only models where the curvature (of either sign) of the inflaton potential 
is not very suppressed with respect to H can give rise to a blue |L6| (red ||I7| ) tilted 
spectrum. 

Note that COBE requires the product \km^\ to be extremely small. If we take 
k ~ 10 -10 , motivated by axion physics as discussed above, then this implies in the 
electronvolt range. The requirement of such a small mass leads to several interesting 
requirements on the model. We envisage that at the Planck scale the <fr mass is equal 
to zero. This can be naturally accomplished within the framework of supergravity 



no-scale models [ 1S[| , where some (not necessarily all) of the SUSY soft masses are 
predicted to vanish, but with non-zero and universal trilinear coupling parameters. 
However the high energy value of will be subject to radiative corrections which 
are very small, being controlled by the small coupling k. In our model the radiative 
corrections at <fi ^> <p c arise from loops of the scalar and pseudoscalar components 
of the complex iV field, which are split by soft SUSY breaking terms, and by their 



fermionic partners. The result may be easily obtained to be ^(J : 

where \x is the modified minimal subtraction scale and we have assumed <fi ^> <p c . If 
we take fi 2 = 4k 2 <f) 2 ~ ITeV 2 to remove the logarithm then use the renormalisation 
group (RG) to run the mass from the Planck scale (where it is zero) down to this 
scale then the radiative corrections result in a correction to the mass of: 

S=i-lW, (18) 

with c a constant of order 1. The COBE constraint, km^ ~ 10 -18 , together with the 
naturalness requirement that the radiative correction is of order the mass itself, will 
now translate into, 

leading to a mass in the eV range. Notice that RG equations result in a negative 
mass squared which would, in the absence of any other correction, appear to favour 
the inverted hybrid inflation scenario. 

However there is another potentially large contribution to the <f> mass coming from 
the vacuum energy V(0) which breaks supersymmetry and which drives inflation. On 
general grounds, whenever local SUSY is broken by non-zero F-terms in the visible or 
hidden sector, the Kahler potential will generate soft SUSY breaking parameters in 
the observable sector which leads to the so-called 77-problem. All the scalar fields in 
the observable sector will pick up masses of order of the Hubble constant [|| ^T], f22 



where H 2 « Vo/Mj>, due to the presence of an exponential factor for the Kahler 
potential in front of the potential for the observable sector. The inflaton, unless 



fine-tuning or specific forms of the Kahler potential are assumed |[23|| , will also get 
this type of mass, making it difficult to satisfy the slow-rolling constraint 77 1 - 
A linear superpotential for the inflaton is a special case that can avoid the problem 



H even for minimal Kahler potential, but requires an input mass parameter. One 
simple way out of the problem is to consider V(0) not having an F-term origin, 
but instead originating from some D-term \ ! ZM. Most F-type breaking supergravity 



models assume than the D-term contribution to the potential vanishes when the fields 
in the hidden sector get their VEVs. This may not be necessarily the case if these 
fields are charged under some Gh hidden or visible gauge symmetry, such as a gauged 
(possibly anomalous) U(l) symmetry. However the identification of V(0) with some 
D-term seems problematic, at least within the framework of string theories, since the 
mass scale of the parameter £ predicted by string theories is too large compared to 
that required by hybrid inflation [fffij , and in the present model this problem is made 
worse due to the smallness of V(0). 

The problem of the origin of V(0) in our model is in fact the cosmological constant 
problem and the same problem besets the MSSM where the explicit potential at the 
global minimum does not vanish, Vmssm(< Hi >, < H 2 >) ^ 0. In order to obtain 
a small cosmological constant in the MSSM one is forced to add by hand a vacuum 
energy V(0) to accurately cancel it as in Eq.||. Here, as in the MSSM, we do not specify 
the origin of V(0), but simply add V(0) by hand. The solution to the cosmological 
constant problem almost certainly lies beyond supergravity, and probably beyond 
string theory. What is crucial for the success of our model is that during inflation all 
the F-terms vanish, and this condition is in fact satisfied by all the explicit terms in 
our model. All that we require of V(0) is that it does not originate from the F-term of 
a supergravity model. Such a vacuum energy which has the characteristics assumed 
here, namely that it does not originate from the F-term of a supergravity model, and 
has at least the possibility of solving the cosmological constant problem has been 
proposed within the context of quantum cosmololgy [p6fl . 

Having fixed the value of the Yukawa coupling, the vacuum expectation values are 
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then around the scale ~ iV ~ 10 13 GeV. More specifically, from Eq. (|T9| ) we have, 

1/3 



10 13 (l0 10 k) GeV. (20) 



The value of C slightly exceeds the usually quoted upper bound for the axion VEV, 
derived requiring the cosmic density of axions not to exceed the critical density of the 



universe [27]. However a value f a ~ 10 13 GeV is not strictly excluded when several 



uncertainties entering the derivation of this bound are allowed fl28fl . Recently, it has 
been also argued [29| that a value of f a bigger than 10 12 GeV can be allowed in models 



where the reheating temperature goes below 1 GeV, that is, below the temperature 
at which the axion field begins to oscillate. The point is that during inflation the 
PQ symmetry is broken and the axion field is displaced at some arbitrary angle, 
and it relaxes to zero only after reheating and only below the QCD phase transition 
when its potential is tilted. At this point the dangerous energy stored in the axion 
field is released, but if the reheat temperature is of order 1 GeV then the resulting 
axion density from the displaced axion field will be diluted by the entropy release 
produced by the inflaton decay. On the other hand the inflaton itself may decay 



directly into axions, and this branching fraction must be sufficiently small so that 
the resulting number density of axions at the time of nucleosynthesis is only a small 
fraction of a neutrino species, as we now discuss. 

Immediately after inflation ends the masses of the singlet scalars which correspond 
to the oscillating mode of the fields N and <fi just after they receive their global 
VEVs will be ~ 0(1 TeV). The resulting "inflaton" mode can decay into lighter 
particles, with a width proportional to the coupling k 2 /(4Tr) ~ (1 TeV) 2 / f 2 . As 
discussed in Ref. [|29| in order to avoid conflict with nucleosynthesis in models where 
the reheat temperature is of order 1 GeV one requires the branching fraction of the 
decaying inflaton into axions not to exceed about 10%. As pointed out |29| the 
inflaton coupling to stops m 2 / f a may dominate over the coupling to axions M| / f a if 
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m| > M|, providing that the stop mixing results in a kinematically accessible light 
stop mass eigenstate mj~ < M|. In this case the inflaton may decay predominantly 
into stops. The resulting decay rate may be estimated as: 

-^M ~ l(T 8 eV% (21) 

which is quite suppressed with respect to, 

H~Y^^^~lMeV, (22) 
3M P 3M P ' v ; 

The reheating temperature is given by pi: 



TnH-OMg^JrtMp, (23) 



where g* is the number of effective degrees of freedom at reheating, and is the 
width of the inflaton decay. Conversion of the vacuum energy to thermal radiation 
through the decay of the inflaton <fi into light particles will be quite inefficient, because 
r^, <C H. This gives a reheating temperature Trh ~ 0(1 — 10) GeV. 

Despite its low value, the reheat temperature is high enough to preserve the stan- 
dard scenario for nucleosynthesis, Trh > 6 MeV, although quite far to allow elec- 
troweak baryogenesis. Moreover, any pre-existing baryon asymmetry is likely to be 



diluted during inflation. Nevertheless, as has been pointed out |29|, |17|| , the amount 
of baryon asymmetry needed might be produced directly by the decays of the infla- 
ton. For this mechanism to work we require the presence of baryon-number violating 
operator in the superpotential, type X'^U^DjD^. As discussed the inflaton can de- 
cay predominantly into light stop squarks, and the subsequent decay of the stops 
into two down-type quarks from this R-parity baryon number violating operator will 
generate baryon-antibaryon asymmetry. Other mechanisms, like Affleck-Dine type 
baryogenesis J3l|, might also work. 



We now turn to the question of origin of the extremely small values of the couplings 
A and k which are required in this model, and to whether this might be understood 
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in terms of a deeper theory such as string theory. In fact it has been argued |32|, ^3 
that an approximate PQ symmetry that solves the strong CP problem can arise from 
superstring models with exact discrete symmetries after compactification. Let us 
assume a Z 3 x Z 5 discrete symmetry, resulting from some string compactification and 
introduce singlets M, M, with the fields transforming as in Table 1. All the fields 
are supposed to originate from some 27 and 27 representations of Eg apart from 
which is taken to be a singlet of E & . To be precise one can assume that Hi, H 2 , N, M 
originate from 27's while M originates from a 27. The superpotential is given by, 

,>tttt,t M M ,,, Ar2 M 2 (MM) 3 (NMf(MM) 2 

where all the Yukawa couplings can be assumed to be of order unity, and we have 
included only the leading physically relevant terms. Note that terms such as N 3 are 
forbidden by the E e gauge symmetry since the product of three 27's does not contain 
three singlets. Here the E§ gauge symmetry is assumed to be broken at the string level 
by for example Wilson line breaking [B3J, so questions such as doublet-triplet splitting 



are addressed at the string level. We assume that the M, M fields radiatively generate 
VEVs, as a result of a radiative mechanism due to their soft masses m 2 becoming 
negative, stabilised by F-terms arising from the above operators with resulting VEVs 



/ mMjA 
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(M) = (M) ~ [~^J ~ W GeV (25) 

As a result of these VEVs we recover the two terms of the superpotential given in 
Eq. ©, with 

A ~ A — 5- ~ A 10" 10 , 26 

M 2 

k ~ k' — T ~ fc'10" 10 , (27) 
MJ. v ' 

The U(1)pq symmetry is explicitly broken by the higher order term proportional to 
d. As discussed |32) such higher order terms contribute an explicit axion mass Ami 
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//. U, X o M U 

3 



Zq a 2 a 2 1 a a a 



1 1 (3 p 3 P 2 



Table 1: Z 3 x Z 5 charges for the chiral supermultiplets. 
which tilts the axion potential slightly, and perturbs the 9 angle by an amount 



Am 2 f 2 

A6> ~ . (21 



In order to preserve the PQ solution to the strong CP problem we require A8 < 10 -8 . 
Setting d — 1 and including a trilinear coupling A = 10 3 GeV, the above operator 
leads to A9 ~ 1CT 11 thereby preserving the PQ solution to the strong CP problem in 
this model. Note that the question of domain walls, both coming from the breaking of 
the discrete string symmetries, and those associated with axion domain walls, which 



was discussed in detail in the second reference in [p2]1 , does not arise in our model 
since they are both inflated away. 

To summarise, we have seen that a simple variant of the NMSSM, called 0NMSSM, 
involving two singlets N, <p but the same number of parameters as in the original 
NMSSM, opens up the posibility of solving the strong CP problem and the fi prob- 
lem, as well as providing also a mechanism for hybrid or inverted hybrid inflation 
in the early Universe, neatly side-stepping all domain wall problems. The smallness 
of the couplings A and k can be understood in terms of a string motivated discrete 
symmetry, while the smallness of the mass implies a no-scale supergravity origin 
for this parameter. We do not specify the origin of the vacuum energy V(0) which is 
necessary to drive inflation, and lead to an acceptably small cosmological constant. 
However we do require that the source of the vacuum energy not be the F-term of a 
supergravity model which would lead to an unnacceptably large mass. The mag- 
nitude of the VEVs are of the correct order of magnitude for axion cosmology, albeit 
on the upper edge of the allowed range. However the model may provide its own cure 
since it has a low reheat temperature of around 1 GeV, and the entropy produced by 
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the inflaton decay partially dilutes the cosmic axion density. Nevertheless, with f a 
just on the upper corner, one expects axion dark matter in this model. In conclusion 
the 0NMSSM has many interesting features and solves several outstanding problems 
of particle physics and cosmology. Like N, the inflaton <ft resides in the visible sec- 
tor of the theory, and will mix with the two Higgs doublets leading to the exciting 



possibility of experimentally observable effects [35 
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Note added 

Since our manuscript was resubmitted two further points have been brought to 
our attention. The first point raised by D.Lyth ( |hep-ph/9710347D is that M.Dine 



( |hcp-th/9207045|) has showed that the model of Casas and Ross is not viable due to 



the presence of additional soft operators which are allowed by the discrete symmetry 
and which violate the PQ symmetry at an unacceptable level. In our model the 
dangerous Dine operators such as N M* MN* , and similar higher order operators, are 
not present due to the absence of the N field, and so our model is exempt from this 
criticism. The second point emphasised by A. Riotto (private communication) is that 
in a certain class of no-scale supergravity models Jl8| (those in which a Heisenberg 
symmetry is present) the inflaton receives no mass of order the Hubble constant 



thereby solving the r] problem p3[ . Since we already invoke no-scale supergravity to 
explain the masslessness of the inflaton, it is natural to appeal to this mechanism 
in our model. Radiative corrections to the inflaton mass during inflation would be 
negligible due to the smallness of the couplings k and A. This of course opens up the 
possibility that the vacuum energy during inflation originates from F-terms after all. 
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